A potential sterile neutrino search using a two-reactor/one-detector configuration 
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There is an observed deficit of about 6% in the expected rate of anti-neutrino interactions when 
averaging over many diflferent reactor experiments. While the significance of the deficit is low (98.6 
% CL), there is speculation that a non- interacting "sterile" neutrino could be the cause. In this 
paper we explore the possibility of a two-reactor/one-detector experiment at intermediate distances 
(100-500 meters) to look for a sterile neutrino in the mass range implied by this deficit. We also 
investigate the potential sensitivity of an existing reactor experiment (Double Chooz) which has a 
single Near Detector at distances of 351 m and 465 m from two reactors of identical design. We 
conclude that Double Chooz could investigate sterile neutrino in the Am^ range of 0.002 to 0.5 eV^ 
over 5 years of near detector running. 



I. INTRODUCTION 

There are global hints of the possible existence of ster- 
ile neutrinos from accelerator experiments [1-4], reactor 
experiments [5] and cosmological measurements [6, 7]. 
While the existence of one or more sterile neutrinos is 
not the only possible explanation for these results, it nev- 
ertheless becomes interesting to devise new experiments 
that explore this possibility. In this paper we investigate 
the sensitivity of an existing experiment (Double Chooz) 
to detect the oscillations of electron flavor neutrinos into 
sterile neutrinos, and propose a new experiment opti- 
mized for this search. 

There have been nearly twenty reactor neutrino exper- 
iments at distance of 10-100 meters from reactor cores. 
In a recent review paper [5] the expected rates were 
re-calculated using up-to-date reactor anti-neutrino flux 
predictions and a 3-flavor neutrino oscillation hypothesis. 
The authors found that on average there is a ^6% deficit 
in the rate of observed anti-neutrino interactions mea- 
sured to that expected. This has come to be called the 
Reactor Anti- neutrino Anomaly (RAA), and has been 
taken as an indication that the 3-flavor oscillation hy- 
pothesis may not be complete. To follow up this hy- 
pothesis, the authors performed global fits using rate in- 
formation from these experiments assuming a 3 active 
fiavor neutrino and 1 sterile neutrino oscillation model 
(also known as a 3-1-1 model). Folding the results of these 
fits in with spectral shape constraints from the Bugey-3 
experiment [8] , they calculated best fit oscillation param- 
eter values of Am^4 = 1.5 eV^ and sin^(26'i4) = 0.14 [5]. 

In this paper we present a new method of gaining in- 
creased sensitivity to 3-1-1 oscillation parameters. This 
method utilizes two reactors and one detector, referred 
to as the "two-reactor/one-detector" configuration, and 
relies on the differences in the measured spectrum from 
running each reactor singly. Thus, this method is effec- 
tive only in the case of two reactor cores, each of which 
runs a significant period of time while the other core is 
down. This is indeed the case for the Chooz B reactor 
configuration and the Double Chooz Near Detector, now 
under construction. The Near Detector will be 351 m 



and 465 m from two 4.25 GWth power reactors. With 
the analysis technique described here, this configuration 
will lead to sensitivity in a region of Amf4 and siv?{29u) 
that has not been previously explored. 



II. SHAPE ANALYSIS OF STERILE 
NEUTRINOS 

To help give clarity to the discussion, a simple example 
is described using the two-neutrino oscillation formalism. 
The two-neutrino survival probability can be described 
by the following formula: 



Pee = 1 - sir? {26) sir? 
= sin^ (/3L) 



Am^L 



(1) 



where a = sin(26') and l3{Ej^^) = Am? / AEj^^. For a single 
detector measuring the disappearance of v^s from a sin- 
gle reactor, the sensitivity to Am^ and sit? {20) depends 
on the reactor-detector distance, L. For an appropriate, 
L, the visible Ey^ spectrum will be distorted due to the 
L/Ep^ dependence of the survival probability. It is possi- 
ble to visually enhance this spectral distortion by forming 
a ratio using an independent Ej^^ spectrum obtained at a 
different distance where sensitivity to Am^ and sin^(20) 
is also expected. In this paper we present a mathemat- 
ical formalism for for a two-reactor/one-detector config- 
uration. It is expected that many (but not all) detector 
systematics will cancel if the single-reactor running peri- 
ods are roughly equally interspersed. This is discussed in 
more detail in section III . In addition, although the de- 
tector backgrounds are not expected to vary with reactor 
running period, we assume for the formalism that the rel- 
evant backgrounds have been subtracted. We note that 
this work was initially motivated to explain the tension 
between the rate measurements and spectral measure- 
ments found in the results of the Double Chooz, Daya 
Bay and Reno experiments in fitting for , which is yet 
not fully understood[9-12]. 

There are four distances (three of which are inde- 
pendent) that can be defined for the two-reactor/one- 
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detector configuration: 

(a) Li = distance from detector to reactor Ri 

(b) L2 s distance from detector to reactor R2 

(c) L2-1 = L2- Li 

(d) Li+2 = Li+ L2 

Using the form of equation 1, the ratio of the energy 
spectra obtained from the reactors i?i and R2 can be 
written as: 



^ l^a^sin^(/3£i) 

l-a^sm'{PL2)' ^' 

This ratio can be expressed in a more useful form by 
multiplying the numerator and denominator by a factor 
of 1 + sin^(/3i2)- After simplifying, the expression in 
equation 2 becomes: 



ee 

pR2 ~ 

l + a^ sin (^^2-1) sin (/3L1+2) - sin^(/3ii) sin^(^L2) 
1 - a4 sin'*(/3i2) 

(3) 

where the spectral distortions due to the mixtures of 
Li and L2 in sine functions are clearly revealed. The 
sin (/3i2-i) sin (/3ii+2) creates modulations on the 
spectral shape at both higher and lower frequency than 
the modulation produced only by baseline Li or ^2- In 
the limit that a is small, equation 3 may be expressed 
as: 



pRi 

ee ^ 

l+[l~a^sm^{l3L2)\ [a^ sin (/JLj-i) sin (/SLi+s)] ^ ' 
+ 0(6) + ... 

where the ratio is a function of the survival probabil- 
ity of antineutrino at a baseline of L2 convoluted with 
modulations of the spectral shape at baselines Li_2 and 
Li_|_2. In L/E space these represent different frequencies. 
Since reactor experiments have most of their events in 
the neutrino energy range of roughly 2-9 MeV each fre- 
quency will have a sensitive "sweet spot" in p. Of course, 
complication comes in with the convolution of the detec- 
tor energy resolution function with the measured visible 
energy, which will reduce contrast in the high frequency 
features in the ratio spectrum. These detector effects will 
be covered in the next section, where we investigate the 
sensitivity of this technique using the Double Chooz near 
detector. 



III. DOUBLE CHOOZ NEAR DETECTOR 
SENSITIVITY 

The Double Chooz (DC) near detector is now being 
constructed at Li = 351 m and L2 = 465 m from Chooz 
B reactors Ri and i?2- A five year run time for this de- 
tector (assuming a down cycle of 15% per reactor, which 
was the case for their most recent publication [10]) will 
yield 548 days of data taken with only either Ri or R2 
on. Each reactor operates at 4.25 GWth, giving expected 
rates of about 230 v^s per day from Ri and 130 VeS per 
day from R2. 

For comparison sake, the DC far detector has been 
currently taking data for 2 years [9, 10]. The distance 
between this detector and reactors (near, far) makes the 
technique developed in the previous section statistically 
limited due to the 1 /i^ fall off of neutrino intensity. The 
measured rates are about 28 V^s per day from i?i and 
23 VeS, per day from R2. 

The DC near detector can offer sensitivity to sterile 
neutrinos not only by making a shape measurement at 
Li and L2, but also by identifying the interference terms. 
This is partly due to the fact that there is a frequency 
modulation at a distance of L2-1 = 114 m, a region which 
has already been probed by previous reactor experiments 
located at a distance of L 100 m from the source. It 
is also important to note that the strength of this ratio 
technique with the DC near detector will not come from 
the sensitivity to the disappearance rate, but from dis- 
crimination of the neutrino oscillation spectral shape for 
values of Ajti^ less than 10~^ eV^. This is a region that 
was not accessible to to Bugey-3 spectral shape analysis. 

The formalism developed in section II is applicable 
when there are no backgrounds present. In actual exper- 
iments, backgrounds such as ^Li and fast neutrons will 
have small contributions which must be taken into ac- 
count. By performing a two-reactor/one-detcctor study 
of the backgrounds, complimented by the availability of 
reactor off data when there are only two cores, we as- 
sume a good understanding of these quantities can be 
obtained. In turn, these backgrounds may then be sub- 
tracted from the far and near reactor data with a small 
loss of sensitivity compared to statistical and detector 
resolution uncertainties, as the background subtraction 
is near fully-correlated since it is simply the same detec- 
tor but a different running period. 

Figure 1 shows the calculated visible energy spectra of 
the i?i and R2 reactor data for values of {l^m\^=Q.b eV^) 
without energy smearing. In the ratio analysis of the 
Ri and R2 data, the following systematic uncertainties 



TABLE I. Distance from Reactor Ri and R2 to the center of 
the detector from geodesic surveys. 





L2-1 


Li 


L2 


Ll+2 


Near (m) 


114 


351 


465 


816 


Far (m) 


117 


998 


1115 


2113 



3 



Expected spectra after applying oscillation and core evolution 

> 2000^ 




J 10 



Energy (MeV) 

FIG. 1. Near (-Ri) and far (-R2) reactor data for the DC near 
detector before smearing of energy resolution witli Am^^^ = 
0.5 eV^ and sm'^{26r,ew) = 0.12. 



on the observed spectral shapes have been taken into 
account: 



• Energy resolution: (7 ± l)%/^E[MeV] from DC 
far detector). 

• Detector stability: estimated at 1% from the DC 
far detector [9] measurements of the stability of the 
Gd capture gamma peak. 

• Reactor core size: 3.47 m diameter for the Chooz 
B reactors. We have assumed the neutrinos start 
randomly inside this core, which is a worst case 
assumption. 

• Fuel loading: contributes < 0.01 % uncertainty, 
based on studies where we put in extremes of fuel 
loading uncertainties. 

The combined effect of these uncertainties to the sensi- 
tivity of sterile neutrino oscillations are shown in Figure 
2. In fact, they have very little impact on the exclusion 
domain of this technique with the DC near detector after 
5 years of detector operation. The analysis is is clearly 
dominated by statistical uncertainty. The spectra from 
the re-evaluation of the Bugey-3 shape discrimination [5] 
is also included as a comparison. 



IV. DISCUSION 

It should be noted that the mathematical formalism 
developed in section II is generic and can be applied 
to any two-reactor/one-detector experiment. The pos- 
itive discovery of sterile neutrinos will require strong ev- 
idence and as such a strong shape discrimination and 
good statistics are paramount to any new proposed ex- 
periments. 

There are a few guidelines required to optimize this 
technique for new experiments. First of all, the distances 
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FIG. 2. Exclusion domain of Arn„e^ for the case of the 
DC near detector. Sensitivity to sterile neutrino oscillations 
via spectral shape discrimination disappears for values of 
Am^g„ less than 0.008 eV^ and greater than 0.16 eV^ at 
sin2(26l„e») = 0.12 at 90% CL. 



between the detector and the far and near reactors cannot 
be equal since the interference term would vanish. When 
taking the ratio of two spectra, the strength of shape 
discrimination is correlated with the difference of baseline 
{L2-1) and the addition of the baseline {L2+i). One can 
envision an experiment with L2-1 = 10 ~ 15 m that can 
devised such that the ILL region can be probed, however, 
the reactor core size will add greater uncertainty at short 
distances and will have to be modeled carefully. 

New experiments must also seek to improve spectral 
shape sensitivity. The principle culprit for the loss of 
shape discrimination power is the energy resolution of a 
detector. Ratio analyses for different values of lAm{^ are 
shown with no smearing applied in Figure 3 and with 
smearing applied in Figure 4. If sensitivity to larger 
Arrii^ is desired, detectors with better energy resolution 
is a vital requirement. 

Currently, there are few experiments which have sen- 
sitivity to Am^gjy below the limits set by Bugey-3. The 
ICARUS experiment can claim sensitivity in this region, 
however it is a i'^ — > appearance experiment and will 
not address the probability of D^. disappearance. If a so- 
lution to the RAA exists in the form of sterile neutrino 
oscillations, such a solution could be probed by reactor D^. 
experiments which cleverly address the many unknowns 
present in reactor flux predictions. In addition, some 
cosmological studies (e.g. those using WMAP-9 data) are 
more consistent with 4 rather than 3 flavors of neutrinos 
{Neff = 3.84±0.40) [7]. When combined with BAG and 
Hq data and other CMB measurements, this same study 
gives an upper limit on the total mass of neutrinos of 
mtot < 0.44 eV at 95% CL. If all of this mass were due 
to a sterile flavor, then the region below (0.44)^ ~ 0.19 
eV^ becomes even more interesting to explore with this 
new technique. 
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FIG. 3. Ratio analysis for value of Arrii^ respectively of 2.0, 
0.5, 0.1 and 0.01 eV^ with no detector smearing. Larger values 
of Ami4 lead to higher frequency of oscillation. 



FIG. 4. Effect of detector resolution on the ratio analysis for 
value of Ami4 respectively of 2.0, 0.5, 0.1 and 0.01 eV'^ with 
a 7%/-\/E[MeV] detector smearing. High frequency features 
are washed out by the smearing of the detector resolution. 
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V. CONCLUSION 

The technique presented in this paper provides a for- 
mahsm that could be used to explore possible sterile neu- 
trino explanation tension that the reactor anti-neutrino 
experiment such as Double Chooz, RENO and Daya Bay 
are observing. A positive discovery would imply that 
sterile neutrinos have a mass less than 0.2 cV^, which 
is consistent with current cosmological results. It is ob- 



served that the distances between the Chooz B reactors 
and the DC near detector will allow sensitivity to a new 
region of Amf^ from 0.002 to 0.5 cV^. The near DC 
detector will provide high statistics and will have the 
power to discriminate spectral features due to active fla- 
vors from those due to sterile neutrinos. Future exper- 
iments could seek to optimize the distances involved in 
a two-reactor/onc-detcctor analysis to increase spectral 
shape discrimination in the Am^ region of interest. 
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